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Abstract 



o 

o 

^ I A semiclassical calculation of particle production by a scalar field in a potential 
is performed. We focus on the particular case of production of fermions by a Nambu- 

<^ • Goldstone boson 9. We have derived a (non)local equation of motion for the ^-field with 

\0 '. the backreaction of the produced particles taken into account. The equation is solved 

^ I in some special cases, namely for purely Nambu-Goldstone bosons and for the tilted 

O i potential U{9) oc rn^O"^. Enhanced production of bosons due to parametric resonance is 
investigated; we argue that the resonance probably disappears when the expansion of the 

. universe is included. Application of our work on particle production to reheating and an 

^ ! idea for baryogenesis in inflation are mentioned. 
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I. Introduction 



Nambu-Goldstone bosons (NGBs) are ubiquitous in particle physics: they arise 
whenever a symmetry is spontaneously broken. If there is additional explicit symme- 
try breaking, these particles become pseudo Nambu-Goldstone bosons (PNGBs). In this 
paper, we consider particle production by Nambu-Goldstone bosons as they rotate about 
the bottom of the 'Mexican hat' potential, whether or not there is a tilt around the 
bottom (i.e. with or without explicit symmetry breaking). 

The Nambu-Goldstone bosons, hereafter called 9, are assumed to couple to fermions; 
thus as the 9 field moves it is capable of producing these fermions. Here we perform a 
semiclassical calcuation. The 9 field is treated classically while the particles produced 
are quantized. The backreaction of quantum fermions on the evolution of the 9 field is 
calculated. Our motivation is to lay out a general way to perform such a calculation 
and to carry it out for a specific coupling. The primary results of our calculations in fiat 
spacetime are Eqs. (2.21) and (2.22). We then demonstrate examples and generalize to 
curved spacetime with massive fermions. 

Particle production by Nambu Goldstone fields may have several applications. One 
is the QCD axion. Another is particle production by the inflation field [1] in Natural 
Infiation [2]. Particle production is, of course, important for estimates of reheating in 
infiation. Many models of infiation involve 'slowly rolling' fields that evolve down a 
potential. Subsequent to the 'slowly-rolling' epoch, there must be an epoch of reheating, 
where vacuum energy is converted into the production of radiation energy. The equation 
of motion for the infiaton field is taken to be 

■ ■ dU 
9 + 3H9 + T9 = -—. (1.1) 

d9 

The term r9 is assumed to describe the reheating. P is taken to be the decay rate of the 
infiaton field. This heuristic term really describes much more complicated physics. In 
fact one should accurately calculate the production of particles and its back reaction on 
the infiaton field as it rolls down the potential. Previous work on this subject includes 
References [3,4]. 

In addition, there may be a mechanism for baryogenesis during Natural Infiation. If 
the equivalent of the Peccei Quinn field can be made to carry baryon number, one may 
be able to do baryogenesis as the infiaton is rolling down its potential. This has several 
nice features: i) the same field would be responsible for both infiation and baryogenesis 
and ii) the infiaton could reheat to very low temperatures, perhaps as low as nucleosyn- 
thesis temperatures of ~ MeV. (In fact, for this mechanism to work one would have to 
reheat to below the electroweak temperature to avoid sphaleron destruction, if it is op- 
erative; alternatively the infiaton could generate nonvanishing {B — L)-asymmetry which 
is preserved by sphalerons). The approach is similar to proposals of Affleck and Dine 
[5] and Cohen and Kaplan [6]. It is assumed that the inflaton fleld ^» is complex and 
has a nonvanishing baryon number. The corresponding baryon current generated by the 
classical rolling down of the inflaton fleld is essentially equal to the angular momentum of 
the two-dimensional mechanical motion in the plane (i?e$, /m$). Thus, as the fleld rolls 
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in one direction, it preferentially creates baryons over antibaryons, while the opposite is 
true as it rolls in the opposite direction. (We assume that the decays during reheating 
are baryon number conserving). Thus, no CP violation is required of the particle physics; 
instead those regions of the universe in which the inflaton by chance rolls down the po- 
tential in one direction turn out to be baryon dominated, while those that roll down the 
other direction turn out to be antibaryon dominated. Conveniently, each of these regions 
is inflated to be very large, so that it makes sense for our baryon dominated region to be 
large enough to encompass our observable universe. The requirements to create a specific 
particle physics model for this proposal are restrictive and are discussed below. 

In Section II, we discuss the sample Lagrangian we consider, and calculate the parti- 
cle production for this case. First we perform the calculation in the absence of expansion 
of the universe, and for production of massless fermions. In Section III, we apply the 
results to two specific examples: i) a scalar 9 rotating in a potential without explicit 
symmetry breaking, and ii) the same scalar but now oscillating near the minimum of 
a quadratic potential produced by explicit symmetry breaking; for example, such a po- 
tential may give rise to inflation. In Section IV, we include the effects of expansion of 
the universe and fermion masses. We also discuss the possibility of parametric resonance 
[3,7,8] whereby large numbers of scalars might be produced during reheating in inflation. 
We argue that there is probably no resonance if one includes the expansion of the uni- 
verse. In Section V we discuss the possible application to inflation, particularly to the 
model of Natural Inflation. The baryogenesis model mentioned above is discussed in this 
section. In Section VI we conclude. 

II. Particle Production in a Simple Model 

IIA) Lagrangian: We flrst describe a simple model in which we calculate particle produc- 
tion. Consider the fundamental action for a complex scalar fleld $ and two fermions Q 
and L: 

j d''x^[g^''d^^*dy^-V{^*^) + iQ^^d^Q + iLYd^L+{g^QL + h.c.)]. (2.1) 

Note that Q and L can be any fermions, not necessarily quarks and leptons of the standard 
model. For example, they can be heavy fermions; they may be given some of the same 
quantum numbers as particles in the standard model if they couple to ordinary quarks 
and leptons. For this section of the paper, we will take the intrinsic mass of the Q and 
L flelds to be zero, and will include mass effects in later sections. 

This action is invariant under the appropriate U{1) symmetry. For example, in this 
paper, we will take the Lagrangian to be invariant under 

$^e*"$, Q^e*"g, L^L. (2.2a) 

Eq. (2.2a) is the symmetry we will use for the rest of the paper. 

We did, however, want to point out that a very similar analysis would apply to the 
case of global chiral U{1) symmetry in a Lagrangian with Yukawa coupling g^L'^R^. 
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Here subscripts L and R refer to left- and right-handed projections of the fermion fields, 
'4^R,L = (1 i 75)V'/2- This Lagrangian is invariant under 



ipL ^ e'"/^ijL, V'fl^e-'"^^ $^e^"$, (2.26) 

which is the Peccei-Quinn (PQ) symmetry [9] in axion models. Then the current would 
be J'^ = 7/57^7^ V'- Although we do not explicitly analyze this case, it would be very 
similar to the one we do look at. 

We assume the global symmetry is spontaneously broken at the energy scale / in 
the usual way, e.g. via a potential of the form 



F(|$|) = A 



(2.3) 



The resulting scalar field vacuum expectation value (VEV) is ($) = /e*'^'-'/v2. 

Below the scale /, we can neglect the superheavy radial mode of $ {rriradiai = ^^^^f) 
since it is so massive that it is frozen out. The remaining light degree of freedom is the 
angular variable 0, the Goldstone boson of the spontaneously broken U{1) (one can think 
of this as the angle around the bottom of the Mexican hat described by eqn. (2.3)). For 
simplicity of notation, we introduce the dimensionless angular field 9 = (f)/ f. We thus 
study the effective Lagrangian [10] for 9: 

f2 - - 

^ef f = —d^9d^9 + iQi^d^Q + il-ff'd^L + {gfQLe'' + h.c.)-U{9) . (2.5) 

The global symmetry is now realized in the Goldstone mode: C^f / is invariant under 

Q^e'^'Q, L^L, 9^9 + a. (2.6) 

At this stage, 9 is massless because we have not yet explicitly broken the symmetry. 

With a rotation of the form in Eq. (2.6) with a = —9, the Lagrangian can alterna- 
tively be written 

^eff = ^-^d^9d^9 + iQYd^Q + il^d^L + {gJQL + h.c.) + d^9J^ - U{9) , (2.7) 

where the fermion current derives from the U{1) symmetry; here, = Q'JuQ- 

Explicit symmetry breaking: Our subsequent analysis of particle production applies 
whether or not the symmetry is further broken explicitly. Several options exist for ex- 
plicitly breaking the global symmetry and generating a PNGB potential at a mass scale 
~ A. Models include the the schizon models of [Ref. 12]. Another possibility is the 
QCD axion [13]: dynamical chiral symmetry breaking through strongly coupled gauge 
fields. When QCD becomes strong at a scale Aqcd ~ GeV, instanton effects become 
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important. Chiral dynamics induces a fermion condensate, (V'V') ~ the potential 

for the angular PNGB field becomes 

U{e) = A^[l ± cos6'] . (2.8) 

Such a potential is used in the case of Natural Inflation, although at higher mass scales, 
and will be discussed later. 

Equations of Motion: The equations of motion for the Q fields are 

il^d^Q + d^OYQ +9fL^0 (2.9a) 

and 

id^iQr) - d^eQr -9fL = 0. (2.96) 
The combination of these two equations can be written 

= d^iQrQ) = -igfiQL - LQ) . (2.10) 
The equation of motion for 6 is 

+ U'{9)/f = j^d^J^ = J^^^~ ■ (^-^^^ 

The equation of motion for the L field, which we assume does not transform under 
the symmetry, is 

ird^L=-gfQ. (2.12) 
As above, if we add the equation of motion for the L field, we find 

d^ilrL) = -igfiQL - LQ) . (2.13) 

Particle Production: Here we calculate production of Q and L particles by the an- 
gular 6 field as it rotates around the Mexican hat (which may or may not be tilted). For 
the moment we will neglect expansion of the universe. 

For convenience we will define 

Q = Qoe'^ (2.14) 

so that eqn. (2.9) becomes 

rd^,Qo = i9fe-''L. (2.15) 

We will solve perturbatively the Heisenberg equations of motion presented above. We 
will take the free field to satisfy 7^9^(5m = 0, and make a perturbation expansion 
Qo — Qin + gQi- Eq. (2.15) is solved by 

Qq{x) = Qin{x) + igf J d^yGqix, y)L{y)e-''^y^ , (2.16a) 
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where to lowest order we will take L = inside the integral. Here Gq is the retarded 
Green's function for the Q field and satisfies 5^7'*Gq(x, y) = 5^{x — y). Similarly, the 
solution to Eq. (2.12) is 



L{x) = Lin{x) +igf J d yGL{x,y)Q{y) 

^ L,n{x) + igf [ d'^yGLix, y)Qin{y)e''^y^ . (2.166) 



We take the vacuum expectation value (in Heisenberg picture) 

{vH{x) + u\e{x))/f) = 



-7^{d^Jnvac= 4{{QL-LQ)),ac. (2.17) 



j2 {^H'^^)vac — J 

We perform a semi-classical treatment, where the 6 field is treated classically while the 
fermion fields are quantized and via the equations of motion determine the evolution of 
6. To first order in g, the right hand side of Eq. (2.17) is (LinQin) — 0. Working to 
second order in g, we substitute Eq. (2.16a) into Eq. (2.14), and plug that into the right 
hand side of (2.17) to obtain 

{v'eix) + u'{e{x))/f) = 

J d^yL,n{x)GQ{x, y)LUy)e-''^'^e''^-^ - Qin{x)GL{x, y)Q^n{y)e''^'^^e-''^-^ + h.c.) . 

(2.18) 

Here, Gq and Gl are the Green's functions for the Q and L fields. For the case of 
massless fermions, Gl = Gq. 

We now quantize the free fermion fields, 

QUx) = J2j + vldl^e^'-'^] > (2-19) 

where 6| and d^ are annihilation and creation operators at momentum k and spin s for 
particles and antiparticles respectively. The quantization for the free Lin field is similar. 
Now we evaluate the right hand side of (2.17). The details of the calculation are discussed 
in Appendix A. As shown there, we find that 



{L{x)G{x,y)L{y))^4 



d^p f dH 



(27r)7y 2E. 



^[i{p+iy{y-x)]i _ ^ ^ (2.20) 



and a similar expression for the quarks. 

After much algebra (presented in Appendix A), we find 



4^2 roo pQ 

{V^e + U'{e)/f) = — V / dw^^ / dt'sm{2ujt')sm[e{t + t') - e{t)]. (2.21) 

Jo J-oo 
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As shown in Appendix A, after the oj integral is done and with some algebra, this 
equation becomes 

e{t + t')cosM-e'^{t + t')suiM . 

(2.22) 

Here w tends to infinity and we have defined = 9{t + t') — 9{t). 

Eqs. (2.21) and (2.22) are the major results of this section. To reiterate, these 
equations describe the evolution of the 6 field with production of massless fermions taken 
into account in a semiclassical approximation. So far the expansion of the universe has 
not been included. 

III. Examples 

In this Section we will apply Eq. (2.21) to two examples. First we will consider the 
case of U'{9) = 0. This is the case where there is no explicit symmetry breaking. Thus 
the potential looks like an ordinary Mexican hat, in which every point around the bottom 
is equivalent. For example, there is no tilt (no cosine potential). Another situation in 
which this case would be relevant is the case when the rotation proceeds so far up in the 
Mexican hat that the details around the bottom are irrelevant and can be ignored. 

The second example we will consider is one that would be relevant to inflation, 
namely oscillations around the bottom of a tilt in the potential. In this case U'{9) ^ 
and there is explicit breaking of the symmetry. 

Case I: U'{9) = 0; In this first case, we consider a simple Mexican hat with no explicit 
symmetry breaking. The zeroth order solution to Eq. (2.21) would be obtained by setting 
the right hand side, which is proportional to g'^, to zero. The zeroth order solution is 
9 = const. In other words, the field is simply rotating around in the Mexican hat with 
constant angular velocity. We substitute this ansatz 9 = const, back into Eq. (2.21) to 
see what the corrections would be. The t' integral becomes dt'sm{2ujt')sm{9t') = 
^[S{2uj + 9)- S{2u; - 9)]. Thus Eq. (2.21) becomes 

9+ ^9'^sign{9) = 0. (3.1) 

If the field $ has been arranged to carry baryon number, then the baryon number 
is shifted (via a baryon conserving decay) to the fermions. The baryon number of the 
fermions satisfies d^J^^ = p9 from the equations of motion. Thus riB = f^9, where ns 
is the baryon number carried by the fermions. The change in baryon number carried by 
the fermions is thus determined by the change in 9 via Ans = f'^A.9, similar to Ref. [5]. 
Here, all that has happened is that this mechanism transfers the initial baryon number 
Bi = f'^9i into the quarks. Subscripts i refer to initial values. Thus it is an initial value 
problem: to get the right value today one would need exactly the right value of 9i. We 
are not proposing this as a likely explanation for the baryon content of our universe. 

The solution to Eq. (3.1) is ^ = 2Tr9i/ {2tt + gH9i) and 9 = 2ng-'^ \ii{2T: + gH9i) +9i. 
Here 9i is the initial value of 9. Thus, we have checked that in the short time or small 



{V^9 + U'{9)/f 



hm 

ZTT to— »oo 



dt' 



cos2wt' — 1 
t' 
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g limit, $ ~ const was reasonable. We wish, however, to remind the reader that the 
expansion of the universe has not yet been included here. 

Case II: Oscillations around the Minimum of a Potential for 6: Here will consider 
the case where U'{6) ^ 0, i.e. there is a tilt as you go around the bottom of the Mexican 
hat. As a simple example we will consider U{9) = as would be appropriate 

near the bottom of a cosine potential. 

Then Eq. (2.22) becomes 



lim 

27r w— >oo 



dt' 



cos2wt' — 1 
t' 



e{t + t')cosA^ -e^{t + t')smAe 



(3.2) 



Again, the zeroth order solution is obtained by setting the right hand side equal to zero. 
Thus as our ansatz we take 6{t) = 0o(^)cosmi?t where mn is the renormalized mass to 
be defined below. We will assume that 9o{t) varies more slowly with time than do the 
cosine oscillations. We will consider the case of small oscillations around the bottom of 
the potential. Then we can take cosA^ 1 and neglect 9'^ in Eq. (3.2). [Note that the t' 
integral should really start from a finite time at which small oscillations begin; here we 
will approximate the lower limit of the integral as ^^^j^j^/ ~ — oo.] Then the t' integral 
becomes 
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27r2 



lim / 



dt' 



cos2wt' — 1 



t' 



[-mji9o{t)cosmR{t + 1')] = 



-^-^ lim 
27r w-^oo 



m%^9osinmRt 



I 



dt'^—^^l^— + m^9osiamRt J dt' 



/dt' 
— (1 — cos2wt')cosmRt' 



smmRt' 



cos2wt' 



(3.3) 



The first term behaves like a friction term, —g'^mR9/{4iT). The second term is zero after 
one takes the limit w ^ oo. The third term is a mass renormalization term. After 
evaluating the integral in this third term as shown in Appendix B, we find that this 

2 

term is £i^m^9\og{2w/mR). We wish to add this term to the second term on the 
left hand side of Eq. [3.2] so that the sum of these terms gives m^9. Therefore we 

2 

define m\[l + log(2t(;/m_R)] = vr? . Then the original integro-difi^erential equation is 
effectively reduced to 



9 + m'-R9 + = , 



where F = g'^ruR/A.'n. The solution to this equation is 

9{t) = 9ie-'^*/^cos{mRt + 6) 



(3.4) 



(3.5) 



where we introduced an arbitrary phase S which is fixed by initial conditions. 

Equation (3.4) (or (1.1)) describes the damping of the external field oscillations due 
to particle production. It was postulated in many papers where the universe reheating 
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was considered. As we have shown it is indeed correct, but our approach alerts us to 
several issues that should be considered further with regard to the calculation of the 
baryon asymmetry. If the spontaneously broken symmetry is associated with the baryon 
number, the baryon asymmetry generated by the decay of the PNGB field was calculated 
[6] as In^l = r/^l^l which gives 

\AnB\ = Tf\Ae\ (3.6) 

Our first caveat is with regard to energy conservation. The initial energy density of the 
field 9 which creates the baryons is p0{ti) ~ prn^Of. At the end some of this energy 
density has been converted to baryons, with energy density psitf) > nsEs where hb 
is the density of the baryonic charge and Eb ~ m is the characteristic energy of the 
produced fermions. Clearly it must be true that PB{tf) < Pe{ti)- If we were to use Eq. 
[3.6] we would see that this requires F < A6m. From the definition of F we see that 
this is satisfied for small values of coupling constant g as long as A9 is not too small; for 
extremely small values of A^, this relationship can never be satisfied. 

Our second caveat is as follows: in making the identification \nB\ — r/^|0|, one is 
equating an operator equation (2.11) with a vacuum averaged equation (3.4). Indeed 
we started with the operator equation which in our case looks like 9 + m^9 — tib/ P- 
Comparison with Eq. (3.4) gives the identification mentioned above. However, Eq. 
(3.4) is not an operator equation but obtained by the vacuum averaging of the operator 
equation (2.11). As we have seen the average value {ub) is not just Vf'^9 but a more 
complicated expression (3.3). This issue should be looked at further. 

Note that in the case of the spontaneous symmetry breaking without any explicit 
one, when U'{9) = 0, the operator equation of motion reads f'^9 = ub and Aub = f'^^9 
(as was mentioned previously). In this case the characteristic energy of the produced 
fermions is 9 and their energy density is of the second order in 9. This agrees with the 
energy density of the creating field 9. Thus the final fermion energy is indeed consistent 
with the original energy in the 9 field. 

The third caveat is that Eq. (3.3) reduced to Eq. (3.4) only in the approximation 
that the lower limit of integration was taken to be — oo. This approximation is good as 
long as there are many oscillations over the course of the integral. For the case of infiation 
the field oscillates many times during the reheating period, and thus this approximation 
is probably reasonable. 

IV. Further Complications: i) Curved Spacetime (Expansion of the Universe), 
ii) Nonzero Fermion Masses, and iii) Parametric Resonance 

Curved Spacetime: So far all our results have been for fiat spacetime. In order to 
include the effects of the expansion of the universe, we now generalize to curved spacetime: 

lf{d^9){dJ)g^^'' + iQV^r^^Q+ 
iLVfJ^'L + d^9QT^Qg^'' - mqQQ - uilLL + gf{QL + iQ) - U{9) 




(4.1) 
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We will consider Friedmann-Robertson- Walker metrics and work in conformal time, ds^ = 
a'^{dT^ — dx^). Since g^i, — a^ry^jy, where rjfj,^, is the flat spacetime metric, by making 
this conformal transformation in the Lagrangian we can reduce the metric to the flat 
one. With this transformation, we can use the ordinary Minkowski space quantization 
for the fields and the usual Green's functions. To simplify we redefine the fermion fields as 
ip ip/a^/'^. Note that the gamma matrices in curved space-time F'* are now transformed 
to normal Dirac matrices Then the action is 



S= d^x 



—rriQaQQ — rriLaLL + gfa{QL + LQ) — a^U (9) 



(4.2) 



where the summations are now done using rj/^^. With this action the equations of motion 
are 

fd^ia'd^e) + d^{Q-i^Q) + U'{e)a^ = , (4.3) 
id^l^Q + d^dYQ - mqaQ = -gfaL , (4.4) 

and 

id^^^'L - rriLaL = -gfaQ . (4.5) 

For the case tuq = tul = 0, we know the fermion Green's functions (for nonzero masses 
there is the complication that ma{T) is time-dependent). In the massless case we can 
repeat the derivation done in flat spacetime previously and flnd that the semiclassical 
equation of motion for the Goldstone fleld is 

{d^{a'df^d) + U'{d)a''/f) = 
— ^a{T) / duu^ / dT'sm{2ujT')a{T + T')sin[^(T + r') - 9{t)] . (4.6) 

Jo J-oo 

Of course, the lower limit of the r' integral should really be some initial time rather than 
— oo. 

Nonzero Fermion Masses: Here we will consider the modifications to the flat space- 
time case when the fermion masses are nonzero. To the Lagrangian in Eq. (2.7) we add 
terms — mgQQ — m^LL. The equations of motion (2.9) are modified to 

il^d^Q + d^e^^Q - rnqQ + gfL = (4.7a) 

and 

id^^iQin - d^dQ^'' + rnqQ -gfL = 0. (4.76) 

Equations (2.10) and (2.11) are unchanged. Again, we wish to calculate Eq. (2.17), and, 
again, {L{x)GQ{x,y)L{y)) is given by Eq. (A. 2). This time we keep the masses mg and 
rriL nonzero. After we perform the J d^y integral, there is a term inside the remaining 
integrals: {I ■ p — rnqmL) / {l"^ - rriQ) = {EiEp - P - rnqmL) / {Ef - P - uiq). There are 



10 



poles at El = iy/^ +™q- Again, only the Ei < part gives a nonzero contribution. 
Our result is 

v^9 + u'{e)/f = 

"5 / "^'^ /I ^'''^"[^^^ + Ei)t']sm[e{t + t') - 0{t)]^^^^^^^^^^. (4.8) 

Here Ef = P + mg and E^ = + m\. Eq. (4.8) reduces to Eq. (2.21) when tuq = 
rriL — 0. 

For the case of ^ = const (Case I considered above), we can see that only particles 
with masses mi +1712 < 9 can be produced in perturbation theory, according to Eq. (4.8). 
If we do the t' integral in ^ = const case, we get 5{d — Ei — E2). For 9 < mi + m2 this can 
never be satisfied, there is no particle production, and ^ = exactly. For 9 > mi + m2, 
this delta function can be satisfied for some momentum, and particles are produced. 

The question remains what happens if one looks beyond perturbation theory, par- 
ticularly in an expanding universe. In the case of e+e~ production by a slowly varying 
electric field, a nonperturbative contribution to particle production exists for the case 
where the oscillation frequency u is less then the electron mass the result is that the 
production is exponentially suppressed (the effect oc exp[—const{me/Lj)]) but nonzero. 
Although we have not found such contributions here, they may exist (see also ref. [3]). 

Parametric Resonance: Recently Kofman, Linde, and Starobinski [7] (see also the 
work of Shtanov, Traschen, and Brandenberger) [8] noticed that parametric resonance 
may greatly enhance the production of bosons during reheating in infiation; one can 
interpret this as the formation of a Bose condensate. In this paper we have been primar- 
ily considering the production of fermions, for which there is no parametric resonance. 
However, we should also consider the production of 9 bosons themselves by the classical 
6'-field. 

We will consider particle production during the reheating phase of Natural Inflation. 
As our potential for the PNGB field we take U{e) = A^{1 - cos^). Then the mass of 
the PNGB field is = A^//^. For Natural Infiation, / ~ nipi and A ~ 10^^ GeV, so 
that m ~ 10-*^^ GeV. We will neglect coupling to fermions in our study of the possibility 
of resonance, and include only coupling of the field to itself. At first, we will neglect 
expansion of the universe and see that parametric resonance does indeed exist for a few 
particular choices of wavenumber. Then we will include expansion and argue why we 
believe that the resonance disappears. We have not performed a complete analysis of the 
equation in the case of expansion, but for the case of Natural Inflation the arguments 
are quite robust. We suspect that the disappearance of the resonance in an expanding 
universe may happen in other cases as well. We leave investigation of this effect in other 
cases to future work. 

In flat spacetime, the equation of motion for the PNGB is then 9 + U'{9)/ f"^ = 0. 
For our choice of potential, this becomes 9 + m^sin^ = 0. For small oscillations about 
^ = 0, we find the solution to this simple equation to be 

9o{t) = 9iSinmt . (4.9) 
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Here, subscript i refers to initial value of the unperturbed solution. Following the ap- 
proach of references [7] and [8], we will now add fluctuations 9 = Oq + 59. We will keep 
terms to first order in 69. Then U' = A'*sin(0o + = A.'^[sm9ocos{69) + sin(5^)cos^o] ~ 
A^[sin^o + 5^cos^o] in the small angle approximation. To first order in 69, after perform- 
ing a Fourier transform and substracting the zeroth order piece, the equation of motion 
for 59 is then 



69 + 



k + cos6'oTO 



69 = Q. 



(4.10) 



Expanding around = and using Eq. (4.9), we have 



777 

69 + + 9lsin^mt]69 = , 



(4.11) 



We define y — 2mt and write sin^mt = (1 — cos2mt)/2. Eq. (4.11) can then be written 



dy- 



■69 + 69 



4;;^ + 4-16^^ + 16^^^°^^ 



= 0, 



(4.12) 



The standard form of the Mathieu equation is 

d'^z 



dy^ 



+ {A + 2ecosy)z = 0. 



(4.13) 



Our fiat spacetime Eq. (4.12) is of the form of the Mathieu equation with A — 
1/4 — ^1/16 + k"^ /Aw?' and 2e = 9^/16. Since we made the small angle approximation 
earlier on, we've assumed 9i < 1, i.e. e <^ A < 1. The Mathieu equation has been studied 
in great depth. It is known that for A — 1/4 + Ai€ and e -C 1, there is no instability if 
I All > 1 (In our case \Ai\ = 2 for k = 0.) Therefore, for k = 0, there is no instability. 
However, for particular values of nonzero k (values for which A ~ n'^/4 where n is an 
integer) there are indeed regions of resonance with 69 growing exponentially in time. We 
refer the reader to literature on the Mathieu equation to see these regions. 

However, now let us include the expansion of the universe. To simplify we will neglect 
here the interaction with fermions. Without fcrmions it is convenient to work in terms of 
the physical time t with the interval ds"^ = dt^ — o?{t)dr^. The relevant change from the 
nonexpanding case will not only be the additional term ?>H9 in the equation of motion; 
rather the important features are i) the redshifting of the lengthscales of the perturbations 
(k ^ k/a) and ii) the fact that the unperturbed solution 6*0 is now different. The equation 
of motion (again, neglecting the fermion effects) becomes 9 + 3H9 + U'{9)/ f"^ — 0. For 
our choice of potential, this equation becomes 9 + 3i70 + m^sin^ = 0. Again, we take the 
small angle approximation sin6' ~ 9. The unperturbed solution for the matter dominated 
expansion is 

9iSni{mt + 6) 



9o{t) = 



mt 



(4.14) 
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The factor of 1/t in the denominator will prove to be an important feature of the expan- 
sion. We assumed here that the 6'-field started to oscillate when the Hubble parameter, 
H = 2/3t, was close to the value of the mass of the field m, as is usually the case for 
inflation. It fixes the initial value of time, ti ~ 1/m. 

This time we take y = mt. The equation for the fluctuations becomes 



^2 SH d 

:Se H -se + 



dy^ m dy 



+ 1 - 7:0f 



rri^a? 2 * y"^ 



S9 = 0. (4.15) 



To eliminate the 59 term we define 9{t) = 9(t)/a^^'^. During the reheating portion of 
inflation, the universe is matter dominated and we take a oc t^/^ (of course after reheating 
the universe is radiation dominated). With this matter dominated expansion, Eq. (4.15) 
becomes 

P 9f sin^y 



69 + 



dy 



56* = 0. (4.16) 



As before, we use s\v?y = [1 — cos(2y)]/2 to write the equation in the form closest to the 
Mathieu equation. This time it is not exactly the Mathieu equation because of the time 
dependence in the denominator of k'^ /a? and because of the factor of l/y^ that came 
fromEq. (4.14). 

We have numerically integrated Eq. (4.12) without universe expansion and Eq. 
(4.15) with the expansion to see what happens to the resonance. As we expected the 
solutions of eq. (4.12) show the resonance behavior for a particular region of parameters 
while solutions of eq. (4.16) do not resonate. The resonance might be excited if the 
oscillations of 9 began when iif ~ 1/tj -C m (not what usually happens in inflation). We 
did not perform a numerical study of the entire range of parameter space, and in principle 
could have missed the particular choices oik/ a that do resonate. Hence we proceed here 
with a simple analytic discussion. 

One can see from simple analytic arguments that Eq. (4.16) is unlikely to lead to 
resonance. When one includes expansion, there are two effects that reduce the instability. 
First, the redshift of the wavenumber, k k/a, quickly moves any wavenumber that 
happens to be in a resonance band, out of the resonance region. In other words, if at one 
time there is an instability on some lengthscale, shortly afterwards this lengthscale has 
redshifted to a value for which there is no instability. Thus it is difficult to see how there 
could be exponential increase in particle production (numerically we could get factors of 
a few, not of 10^). Second, in the long time limit, y ^ 1, both the first and third terms 
inside the brackets become small. Then the equation is simply a harmonic oscillator 
equation with oscillating rather than unstable solutions. The third term, whose negative 
sign could make it responsible for resonance, becomes unimportant for times t > 9i/2m. 
Since the frequency of oscillations m is usually faster than the frequencies corresponding 
to other relevant timescales, such as the timescale of reheating in inflation, the third 
term quickly becomes unimportant. There may be enhancement in the first oscillation 
or two, but it is probably not very large (as above, this assumes that the oscillations 
began when iJf ~ m, i.e. y ~ 1). Instead the solutions quickly become oscillatory, with 
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amplitudes at most slightly larger than those in the non-expanding case. Again, we have 
not performed a complete analysis of the equation (4.16), but we have argued why we 
believe the resonance effects are not very strong here. 

V. Baryogenesis in Natural Inflation 

The inflationary universe model [1] provides an elegant means of solving several 
cosmological problems, including the horizon problem, the flatness problem, and the 
monopole problem. In addition, quantum fluctuations produced during the inflationary 
epoch may provide the initial conditions required for the formation of structure in the 
universe. During the inflationary epoch, the energy density of the universe is dominated 
by a (nearly constant) vacuum energy term p ~ pvac-, ^md the scale factor R of the universe 
expands super luminally (i.e., > 0). If the time interval of accelerated expansion satisfles 
At > 60R/R, a small causally connected region of the universe grows sufficiently to 
explain the observed homogeneity and isotropy of the universe, to dilute any overdensity 
of magnetic monopoles, and to flatten the spatial hypersurfaces (i.e., O — 1). 

The model of Natural Inflation [2] was proposed to provide a natural explanation of 
the required flatness of the potential in inflation. The flatness is achieved by mimicking 
the axion physics described earlier. The inflaton is a PNGB fleld and two different mass 
scales describe the height and width of the potential. The model has several nice features, 
including the possibility of extra large scale power in the density fluctuations, negligible 
production of gravitational waves, and possible tie-ins to particle physics models under 
consideration. 

Here we want to consider an idea for baryogenesis during Natural Inflation. More 
standard ideas such as a) reheating to above the baryogenesis temperature (e.g. elec- 
troweak) or b) baryon violating decays have already been considered [2]. Instead, here 
we consider a model of baryogenesis in which the baryon number is produced as the 
inflaton is rolling down its potential. In this paper we will merely suggest the idea, and 
leave study of the implementation of the idea for future work. 

For this particular idea to work, the inflaton would have to carry baryon number. 
If the inflaton rolls clockwise down the hump in the Mexican hat, then baryons are 
produced; if the inflaton rolls counterclockwise down the hump, then antibaryons are 
produced. In different regions of the universe, there will be these two different kinds of 
behavior. Any one region will be blown up to become very large by the inflation. Thus 
our observable universe, which lies inside one of these regions, had a flfty /fifty chance of 
being made primarily of baryons/antibaryons. CP violation is not explicitly required in 
the Lagrangian, as the sign of the baryon number is determined by the initial conditions, 
namely the direction of the roll of the field. These ideas are very similar to those of [5,6]. 

The reheating temperature in this scenario could be very low. In particular, if 

Tr cheat < Telectroweak , sphalerous do not crasc any baryon asymmetry generated during 
inflation. Of course we also need to return to the standard evolution of the universe at a 
high enough temperature for nucleosynthesis; i.e. Treheat > Tnudeosynthesis- It may be a 
nice feature to have a very low reheating temperature, as many inflationary models are 
very constrained by the requirement of a high reheat temperature. 
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In order for this to work, the $ field must carry baryon number. The current that is 
exphcitly broken by instantons or by whatever else provides the tilt (the cosine potential) 
cannot be orthogonal to baryon number. In that case the baryon number of the $ field 
will be proportional to the angular momentum as the inflaton rolls down: one direction 
of roll will correspond to baryon production and the other to antibaryon production. 
The baryon current carried by the $ field is J'* = i[^d'^^* — Since ($) = /e*^, 

the baryon number density (ub) = (Jo) = /^^, namely the angular momentum of the 
two-dimensional mechanical motion of the $ field in in the plane (i?e$,/m$). As an 
example, in the Lagrangian in Eq. (2.7), we have considered a symmetry whereby $ and 
Q transform whereas L does not. Thus $ and Q could carry baryon number while L does 
not. Q and L would not be ordinary quarks and leptons; rather they would be hidden 
sector particles that could be made to couple to quarks and leptons in such a way that 
Q carries baryon number while L does not. 

There are many constraints that such a model must satisfy. One must be careful 
about the quantum numbers carried by the various fields: namely SU(3) color, the gauge 
group that became strong at scale A and produced the cosine potential in the first place, 
and baryon number. One must also ensure that the present day violation of baryon 
number predicted for ordinary matter is not in excess of observations. Also, we do not 
want the only decay mode to be to baryonic matter of our universe. Somehow there 
must be nonbaryonic decay modes or decays to baryons that remain in the hidden sector. 
Simultaneously satisfying all these constraints is difficult. However, we have by no means 
exhausted all the possiblities, and leave this investigation to future work should the idea 
prove promising enough. As indicated near the end of Section III above, the calculation 
of the baryon number produced in such a model will be performed in future work. 

VI. Conclusions 

As a Nambu-Goldstone boson 9 moves in a potential, it can produce fermions that it 
couples to. A semiclassical calculation of particle production was performed. The back- 
reaction of quantum fermions on the evolution of a classical 9 was calculated for a specific 
simple model, to provide a general framework within which one can calculate production 
of other particles as well. The primary results of our calculations in fiat spacetime are 
Eqs. (2.21) and (2.22). Generalization to curved spacetime with massive fermions was 
discussed. We argued that enhanced production of bosons due to parametric resonance 
is probably not important here in an expanding universe; a more general investigation 
of the effects of expansion on resonance is warranted in the future. We are especially 
interested in the model of Natural Infiation, in which the infiaton is a pseudo Nambu 
Goldstone boson. It may be possible for the inflaton to create baryon asymmetry at the 
exit from inflation simultaneously with the universe reheating. 
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APPENDIX A: 

CALCULATION OF SEMICLASSICAL EQUATION FOR 9 FIELD 



We wish to calculate the right hand side of Eq. (2.17). Using the quantization in 
Eq. (2.19), we find that 

{L{x)GQ{x,y)L{y)) = 
(27r)3 >' 



EE 



s s' 



s -ip-{x-y) 



(27r)3 Ep 

where the second equality follows since {dp dp) = 5as'5^{p — p'). Using 

. f dH exp[-il ■ (x - y)] 
J (27r)* /^7'^ -niQ + ie 

^p^p = (Pfj,!^ - mL)j2mL, we find 

{L{x)GQ{x,y)L{y))^ 



(Al) 



{A.la) 



f d^P f 
J i^TryJ 



dH 

lEr, 



In the massless limit considered in Section II, the trace becomes Tr[l^p^^^'y'^] = 41 ■ p 
and we have Eq. (2.20). The calculation of {Q{x)GL{x,y)Q{y)) is similar. 

Using Eq. (2.20), we now have terms on the right hand side of Eq. (2.17) such as 

J J (27r)' Epr 

where we have taken / —I compared to previous expressions. We now write — {Ei, I) 
and p^ — {Ep, p) so that l-p — EiEp — l -p. Henceforth we will assume no spatial gradients 
in the 9 field, i.e. 9 = 9{t) only. We wiU now use / d^ye^^^P-'^-y^ = {277^6^ {p- V) to write 



(27r)4 Ep 



Ef - P 



^h.c. (A4) 



We now do J dEi and find a nonzero contribution from the pole at Ei = — Note that 
we take the retarded Green's function, which gives nonzero result only for ty < tx- We 
have 



9 

27r2 



-oo 



(A5) 
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Adding all terms of this form that contribute to the right hand side of Eq. (2.17), 
relabeling Ei as a;, and defining t' = ty — tx, we find 



{V^9 + U'{9)/f) = — V/ dw^^^ dt'sm{2iot')sm[e{t + t') -9{t)]. (A6) 

Jo J-oo 

This is the result quoted in Eq. (2.21). 

We can rewrite this in the form given in Eq. (2.22) if we now perform the u) integral: 

1 ^2 



dujsm2ujt' 



= hm ^ 



( cos2wt' — 1 

2t' 



Next we do the t' integral by parts. We will use the notation A6' = [9(t + 1') — 0{t)\. The 
nonvanishing contribution is 



- lim 

8 w^oo 



f" ,,d^\l{ 



cos2wt' — 1 



= —— lim 

8 u)— >oo 



- lim 

8 w— >oo 





— oo 


-oo 



-'i 



cos2wt' — 1 



dt' 



1 
t' 



cos2wt' — 1 



dt'^ 



sinA6' 



^sinA6' 
dt' 



sinA6' , 



(^.8) 



where surface terms have all vanished. We now perform the derivative on sinA^. Then 
Eq. (A. 6) becomes 



lim 

27r w— ►oo 



dt' 



cos2wt' — 1 
t' 



e{t + t')cosAe -e^{t + 1') slnA^ 



This is the result quoted in Eq. (2.22) 



(A.9) 



APPENDIX B: 
CALCULATION OF MASS RENORMALIZATION TERM 
FOR THE CASE OF SMALL OSCILLATIONS AROUND THE MINIMUM 

To obtain Eq. (3.4), we must calculate the following term: 

171^6 ocosniRt / "[I ~ cos2wt']cos'mRt' 
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m%eoCOsmRt J — 



cosniRt' — cosat' + cosniRt' — cos/3t' 



niROQCOsmRt 



t' 



sm{mR + w)t'smwt' + sint(;t'sin(t(; — mR)t' 



{B.l) 



where a = 2w + and P = 2w — rriR. After doing the integral, we find that this term 
is 



cosmut 



log 



/ 2w 



+ rriR 



rriR 



+ log 



2w — rriR 



rriR 



In the limit — > oo, this becomes 

-niROocosmRt log{2w /niR) , 
a logarithmically divergent term that renormalizes the mass. 



{B.2) 



(S.3) 
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